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We show that ordinary and radiative muon capture impose stringent constraints on sterile neutrino 
properties. In particular, we consider a sterile neutrino with a mass between 40 to 80 MeV that has 
a large mixing with the muon neutrino and decays predominantly into a photon and light neutrinos 
due to a large transition magnetic moment. Such a model was suggested as a possible resolution to 
the puzzle presented by the results of the LSND, KARMEN, and MiniBooNE experiments [l[. We 
find that the scenario with the radiative decay to massless neutrinos is ruled out by measurements 
of the radiative muon capture rates at TRIUMF in the relevant mass range by a factor of a few in 
the squared mixing angle. These constraints are complementary to those imposed by the process 
of electromagnetic upscattering and de-excitation of beam neutrinos inside the neutrino detectors 
induced by a large transition magnetic moment. The latter provide stringent constraints on the 
size of the transitional magnetic moment between muon, electron neutrinos and N. We also show 
that further extension of the model with another massive neutrino in the final state of the radiative 
decay may be used to bypass the constraints derived in this work. 

PACS numbers: 14.60.St,13.15.+g 



I. INTRODUCTION 

Muon capture has long served as a powerful probe 
of the weak interaction and of nuclear structure. The 
closely related process of radiative muon capture (RMC), 
in which muon capture is accompanied by a photon 
emission, is greatly suppressed compared to ordinary 
muon capture (OMC), with a typical branching of 
rRMc/r MC - 10~ 5 . The small rate of RMC allows 
for it to place tight constraints on models that would 
lead to increased rates of RMC. In this note we show a 
new application of the OMC and RMC results as power- 
ful constraints on sterile neutrino properties. We argue 
that in the kinematic range accessible to muon capture, it 
provides the best constraints on prompt radiative decays 
of sterile neutrinos. 

Considerable experimental and theoretical effort has 
gone into the study of neutrino properties with medium 
energy proton beam experiments, such as LSND 2], 
KARMEN @1, and MiniBooNE Q. As it is well-known, 
no definitive conclusion can be made about the pres- 
ence or absence of — > v e neutrino oscillations, and 
claimed signals are difficult to reconcile with other mea- 
sured properties of neutrino masses and mixings. While 
it appears not impossible that very small contributions 
to the fluxes of v e might have been overlooked in the 
neutrino production simulations and while some back- 
grounds may not have been fully appreciated at Mini- 
BooNE 5] , it is nevertheless tempting to interpret these 
anomalies in terms of new physics. One such attempt 
was undertaken recently by Gninenko [l|, @ , who con- 
structed a phenomenology-driven model of sterile neu- 
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trinos with sub-GeV mass, a large transitional magnetic 
moment, and a sizable mixing with the muon- type neu- 
trino. The model is designed as a possible explanation 
of some 6] or all "unusual" neutrino signals [l[ at LSND 
and MiniBooNE. 

Sterile neutrinos, among other very weakly interact- 
ing particles, can be a legitimate search object for differ- 
ent high- intensity/high-luminosity experiments (see, e.g. 
PI)- Once produced, a heavy neutrino N with < m„ 
would typically decay into three leptons, and occasion- 
ally to v^. Reference |l| considers a modification with a 
much enhanced transitional magnetic moment that con- 
siderably shortens the heavy neutrino lifetime. To ex- 
plain all the oscillation measurements while evading con- 
straints from K^2 and 7r M2 decays Ref. [l| requires: (1) 
that the heavy neutrino, N, have a mass in the range 
40 MeV < m N < 80 MeV; (2) that its mixing with the 
muon neutrino have a strength 10~ 3 < | CT^^jv | 2 < 10~ 2 ; 
(3) that it decays primarily to light neutrinos and a pho- 
ton; and (4) that its lifetime is tm < 10~ 9 s _1 . Many 
different signatures in meson physics can be entertained 
with such an object and were described in Ref. We 
notice that a sterile neutrino that mixes with the muon 
neutrino and decays dominantly to lighter neutrinos and 
a photon will greatly increase the RMC rate through the 
diagram in Fig. [TJ Our calculations show that the whole 
mass and mixing range suggested by Gninenko is ruled 
out, if the final state of the radiative A~ decay contains 
massless standard model (SM) neutrinos. On the other 
hand, we find that the model modified by introduction of 
another masive neutrino in the final state of the radiative 
decay of A" can pass the RMC constraints. 

In Sec. HH we examine the constraints that RMC mea- 
surements have on sterile neutrinos that decay radia- 
tively. Section IIIII contains estimates of the sterile neu- 
trino electromagnetic production cross sections in a neu- 
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FIG. 1. Diagram that leads to RMC induced by a radiatively 
decaying sterile neutrino N. First, OMC produces an on-shell 
N followed by the decay N — 5* jv. 



trino beam scattering on nuclei through a large tran- 
sitional magnetic moment. We point out that such a 
process provides an additional source of constraints and 
should introduce important modifications to the analysis 
of Ref . [l[ . We discuss the effects that the production of 
sterile neutrinos through OMC on future fi — >■ e conver- 
sion experiments in Sec. |IV] We conclude in Sec. |V] 
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FIG. 2. Left: the helicity asymmetry Alr as defined in 
Eq. @ as a function of mjv. As expected, Alr — 1 for 
mjv = and Alr — > for rriN — > m^. Right: the photon 
energy spectrum [see Eq. J2|] with a — 1 [see Eq. Q] for 
rriN — 60 MeV (solid curve) and mjv = 80 MeV (dashed 
curve). The spectra do not go smoothly to zero near their re- 
spective maximum energies because Alr < 1 for both of these 
masses. The plots are for iV-producing OMC on hydrogen. 



II. RADIATIVE MUON CAPTURE FROM A 
STERILE NEUTRINO 

For simplicity we shall assume the same type of V — A 
interaction governs the transition between the muon and 
the sterile neutrino. Then, we approximate the RMC 
rate induced by the production and subsequent radiative 
decay of a sterile neutrino, N, by the following formula: 

^ £ = |C/^| 2 r Mc(™iv)B(7V^ 7 ^/( J B 7 ) , (1) 

where Tomc (wijv) is the rate of ordinary muon capture 
with a massive sterile neutrino N in the final state, mod- 
ulo the mixing factor U^n- f (E-y) is t ne energy distri- 
bution of the photons, 
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The parameter a is defined by the angular distribution of 
the photon in the sterile neutrino's center-of-mass frame 
through 
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where 9 is the angle between the photon's direction and 
the spin of N . In the most conservative case of a mini- 
mum of new fields beyond the SM, N and v are Majorana 
neutrinos, a — 0, the decay of N is isotropic and the en- 
ergy distribution of photons in the lab frame is flat, cf. 
Eq. ([2]). If the neutrinos are Dirac, the parameter a may 
be in the range —1 to 1, which introduces anisotropy in 
the rest frame of TV decay, and affects the energy dis- 
tribution. In the case of Dirac SM neutrinos with only 



left-handed couplings, a ~ 1 ||. ^-ljj is the helicity 
asymmetry of the sterile neutrinos produced in OMC, 
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with Nl^r the number of left- (right-) handed sterile neu- 
trinos produced. Due to the assumed V — A structure of 
the current producing the sterile neutrino, A l r = 1 for 
friN — 0, leading to no dilution of the photon anisotropy 
with respect to the sterile neutrino momentum. However, 
Alr decreases as tun increases, partially washing out the 
anisotropy and leading to a flatter energy spectrum. We 
show Alr as a function of tojv and the photon energy 
spectrum for a = 1 and mjv = 60, 80 MeV in Fig. [2] for 
OMC on hydrogen. 

We focus on the case of muon capture on hydrogen 
since the theoretical uncertainties due to nuclear physics 
are lessened compared to capture on complicated nuclei. 
Neglecting second class currents, the nuclear matrix ele- 
ment can be written as 



( n \ J w\P) = U n iP2) 



lip 
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with M np = (m„ + m p ) /2 and q = p2—pi- It is necessary 
to calculate the rates for the cases where the muon and 
proton form a spin singlet and a spin triplet separately. 
Using this matrix element, ^ S qmc ( m A r ) an d Tomc (ttin) 
can be found. We use values of the form factors from 
Ref. @ at q 2 = — 0.88m^ + 0.89m^ and assign an error 

of ±20% to gp. In Fig. [21 we show i?^ ng and i?^r lp where 
each is the ratio of the OMC rate for a neutrino of mass 
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FIG. 3. -R^ ng (solid curves) and -R^ lp (dashed curves) as 
defined in Eq. Each set of two curves is generated by 

varying gp by ±20%. 



mjy to the OMC rate for a massless neutrino: 

rising / \ 
T OMC \ m v = m N) 



rjsmg 



r OMC ( m » = °) 



(6) 



and similarly for i?^ lp . 

E452 at TRIUMF measured the rate of RMC on hy- 
drogen for photons with an energy > E^ ut = 60 MeV 
and found a RMC branching ratio LLC 



RMC 



_B T >60 MeV 



= (2.10 ±0.21) x 1CT 8 , (7) 



where r to t is the total muon decay rate. Using this mea- 
surement, we can set a limit on | U^n | 2 B ( TV — > ^v) by 
requiring that the rate of RMC induced by the radiative 
decay of a sterile neutrino with E 1 > 60 MeV not exceed 
this value, once we know the OMC rate for TV production. 

In Ref. [13] photons that originated outside of the hy- 
drogen target were rejected. Therefore, we must multiply 
our estimate of the RMC rate in Eq. ([1]) by the proba- 
bility, P(mjv), that the decay TV — »• ju occurs in the 
target. To estimate this probability, we calculate the av- 
erage distance traversed by a sterile neutrino, d av , from 
a muon capture event to the edge of the target (a cylin- 
der with a 16 cm diameter and a 15 cm height) assum- 
ing that the captures happen uniformly throughout the 
target and that the sterile neutrinos are emitted isotrop- 
ically with respect to the detector orientation. We find 



13.7 cm. The probability is then estimated as 



P (that) = 1 — exp 



7VT 



This probability is approximately 18% for = 40 MeV 
and 43% for 80 MeV if r (TV -> jv) = 10~ 9 s. 

To accurately compare with experiment, the relative 
fractions of different \ip bound states in the target used 
in Ref. [Io[ must be taken into account. Muons are cap- 
tured in the atomic singlet and triplet states. In the 
target used, liquid hydrogen at 16 K, the atomic triplets 
quickly transition to the atomic singlet. These atomic fip 



singlets tend to form molecular p/ip states that can either 
be in the ortho (proton spins parallel) or para states (pro- 
ton spins antiparallel) . This results in the muons being 
captured from a mixture of atomic singlet and molecular 
ortho and para states. Reference [10( found the relative 
occupancies to be 
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0.061, /para = 0.085, /ortho = 0.854 . (9) 



The molecular OMC rates can be related to the atomic 
ones, 
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where the approximately equal sign in the ortho rate re- 
flects the fact that the ortho molecular state is thought 
to be dominantly spin-1/2 with little spin-3/2 contribu- 
tion [ll[. The wavefunction corrections are 7p = 0.5733 
and 70 = 0.500 QJ]. The total OMC rate is then 
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Knowing the atomic singlet and triplet OMC rates as 
functions of rriN and the relative abundances of different 
muon-containing modifications of hydgrogen, we calcu- 
late Fomc (tun)- This rate, the relationship in Eq. ((T|), 
the measurement in Eq. ([7|1. the energy spectrum in 
Eq. @, the probability in Eq. ©, and the assump- 
tion of a lifetime r (N — > 72/) < 10~ 9 s allow us to limit 
\Ufj, N \ 2 B(N -> ^v). This limit is shown in Fig. U for 
a = —1,0,1. For a branching ratio B (N — > ^v) ~ 1, 
this implies that | C/ m jv | < 7.9 x 10~ 4 over the range 
40 to 80 MeV for a = 1, slightly below the mixing re- 
quired in Ref. fl. For a = 0, -1, with B(N -fyu) ~ 1, 
the limit over this mass range is stronger, iL^jvj 



< 

(3.3,2.1) x 10~ 4 respectively. We note that this is a 
conservative upper limit since we have simply required 
that the number of RMC events due to the radiative 
decay of the sterile neutrino N not exceed the total 
number of events seen. In practice, TV-induced RMC 
occurs along with standard RMC, and requiring that 
the sum of the two contributions agrees with experimen- 
tally observed rates would additionally strengthen the 
bound by a factor of ~ 2 — 3. We also notice that our 
estimate of the OMC rate is on the slightly low side 
(Tomc = 687 — 711 s _1 as tun — > with our range 
(8) of nuclear form factors as compared to the experimental 
value of rQ^| c = 725.0 ± 17.4 s _1 [13| ) which also makes 



our limit on |t^jv| more conservative. 

Our calculations also show how the RMC limits can be 
weakened or avoided at the price of slight modification of 
the orginal model in Ref. 1]. In order to do that, one can 
further "downgrade" the photon spectrum by assuming a 
massive neutrino TV' in the final state of the radiative de- 
cay, TV — > TV'7. The photon energy is then decreased by 
the factor (l — m 2 N , /mff) . If the TV' is massive enough, 
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FIG. 4. The range of limits on ]f/ M jvj 2 B (N -»■ ju) with m„ = 
for a — — 1, 0, 1 (from bottom to top) implied by the RMC 
rate measured in Ref. [13] for a sterile neutrino of mass mjv 
with lifetime r (N ju) < W~ 9 s. The ranges are generated 
by varying gp by ±20% and the measurement from Ref. 
by ±lcr. The a = 1 case is least constrained since photons 
are then preferentially emitted opposite the direction of the 
sterile neutrino's momentum, cf. Eq. ([2]). The limits weaken 
for smaller mjv since the probability of the radiative decay 
occurring within the target in Eq. (JSJ) decreases. 



the maximum photon energy could be below the experi- 
mental cut of 60 MeV. For to^ = 60 MeV, this requires 
to at- > 35 MeV. 

One should also be aware, of course, of the discrepancy 
in the RMC-extracted value for the gp form factor com- 
pared to the chiral perturbation theory calculations. The 
experimental value for gp from the standard RMC seems 
to be in excess of the standard model (SM) prediction at 
a ~ 30% level [H]]. We note in passing that one could 
speculate that such discrepancy could also originate from 
the radiatively decaying sterile neutrino, with a squared 
mixing angle at the level of O(10 -4 ). 



III. PRODUCING N THROUGH THE 
MAGNETIC MOMENT 

The radiative decay rate of the sterile neutrino to SM 
neutrinos is given by 

8 L VaW rnl 



r (n 



(13) 



where fi\ T is the transition magnetic moment between 
N and the i th lighter neutrino. Factors of m e in the 
above formula are due to a conventional normalization 
of magnetic moments on the Bohr magneton fip ■ If one 
identifies v with the three light neutrinos of the SM and 
assumes a common value of the transition magnetic mo- 
ment between N and each of them, a radiative decay 
lifetime of 10 -9 
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A transition magnetic moment this large raises the pos- 
sibility that these sterile neutrinos could be produced 
by neutrino beams electromagnetically upscattering on 
nuclei. A similar production mechanism was consid- 
ered in [15l | . Muon or electron (anti)neutrinos scattering 
electromagnetically on a target of spin-1/2 with charge 
Ze will produce sterile neutrinos through the process 
v ell Z — > N Z with the cross section (ignoring subdom- 
inant scattering on the nuclear magnetic dipole moment) 
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where M is the mass of the target, T is its recoil kinetic 
energy and F (g 2 ) its electric form factor, and E v is the 
energy of the light neutrino. For a scalar target (e.g. 12 C 
nucleus) this cross section becomes 
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These expressions agree with those found in Ref. [l6| in 
the limit that tun — > 0. The energy of the sterile neu- 
trino is Em = E v — T . At LSND, this process would lead 
to a deposit of electromagnetic energy through v 12 C — >• 
N 12 C — > 7^ 12 C. Using a charge radius for 12 C of 
2.46 fm [13], we show a {y 12 C -> N 12 C) as a function 
of the incoming neutrino energy for rajy — 60 MeV and 
t{N -> ju) = 10- 9 , 10- 10 , 10- 11 s in Fig. assuming 
a common transition moment between N and the light 
neutrinos. Using the flux of and v e from pions decay- 
ing in flight in Ref. Q, one can estimate the number of 
sterile neutrinos produced in the LSND target through 
this process, 
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which would mimic v e C — > e~ N (N refers to nitrogen 
here) events after the radiative decay of the sterile neu- 
trino. The fitted number of events of this type was 18.0, 
indicating that a radiative decay lifetime for N of 10 -9 s 
or shorter is ruled out if there is a common transition 
magnetic moment to the three light neutrinos. 

The production mode due to the transition magnetic 
moment and strong constraints on the overall number of 
electron-like events in LSND may be avoided if transi- 
tional magnetic moments are flavor-nonuniversal. One 
could avoid both the problem of exessive N production 
in the neutrino detector and strong constraints on |t/ M jv| 2 
in Sec. [II] using a hierarchy of magnetic moments, 
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FIG. 5. a [y 12 

r(N^'yv) = 10~ u , 10"" lu , 10 _ii s as functions of the in- 
coming neutrino energy. We have assumed a common tran- 
sition magnetic moment between N and the light neutrinos, 
jUt r /Ms = 9.6 x 1(T 9 , 3.0 x 10 -8 , 9.6 x 10~ 8 respectively. 



and still possibly provide an explanation to at least some 
of the experimental anomalies. If ^ r is decreased by two 
orders of magnitude from the value in Eq. (|14|l and /^J? r by 
a factor of at least several so that T (N — > jv^) < 10 5 s _1 
andT(N -> ju e ) < 10 7 s"\ the total rate of v e C -> e" N 
events could be made consistent with the LSND bound. 
If a large transition magnetic moment between N and v T 
is maintained (/xj" r <1 10 _8 //s), 
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s < r (N -> ju T ) < 10" 



(19) 



the sterile neutrinos would then be produced through 
their transition magnetic moment with v e or while 
they would decay promptly to ^v T . The mixing | t/^jv | 2 
could then be reduced to agree with the limit from RMC. 
While it appears possible to explain the strength of the 
MiniBoone signal that way, a detailed study of the LSND 
signal is a lot more complicated. It would involve some 
nuclear modeling to determine whether the — > N up- 
scattering process could additionally dislodge a neutron 
from the nucleus with large enough probability to look 
like the LSND signal without disturbing other measure- 
ments. It is also not immediately clear without further 
investigation if the angular distribution of photons from 
this production mechanism would be compatible with the 
signal. Finally, one should also assess the constraints on 
the transitional magantic moments of atmsopheric r neu- 
trinos imposed by the super-Kamiokande electron-like 
rates. This model's motivation is purely phenomenologi- 
cal and the question of making a sensible UV completion 
of it must also be addressed. 



IV. SENSITIVITY IN fi ->• e CONVERSION 

Searches for sterile neutrinos with a mass less than that 
of the muon that mix with muon neutrinos are particu- 
larly suited to experiments looking for fi — > e conversion 



in the presence of a nucleus due to the large number 
of muon captures. Large Z nuclei are used as targets, 
greatly increasing the muon capture rate. In aluminum, 
for example, the capture rate and the free decay rate are 
roughly equal, and for larger Z the capture rate begins 
to dominate. 

The proposal for the Mu2e experiment at FNAL [l8| 
envisions around 10 17 muon caputures per year, produc- 
ing roughly 
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sterile neutrinos. Using a 48 Ti target, branching ratios 
of the order 10 -16 would be probed. Mixing angles 



\U»n\ > 
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(21) 



could be observed, with e < 1 describing the fraction 
of decays that would pass experimental cuts, nuclear 
physics, and efficiencies that depend on the details of 
the experimental setup and model under consideration, 
e.g. the lifetime of N. For small U^n, of course, standard 
RMC becomes dominant and information about the spec- 
tral shape is necessary to increase sensitivity. The nuclear 
physics in high-Z materials is also less well known than 
in the case of hydrogen considered in Sec. [TTJ Another 
benefit from such experiments is that sterile neutrinos 
with less exotic decay modes, like ve + e~ , could also be 
explored. 



V. CONCLUSIONS 

In this note, we have shown that RMC can strin- 
gently constrain models of new physics. It is particu- 
larly suited for putting constraints on models of sterile 
neutrinos in which they promptly decay with a photon 
and the massless neutrino in the final state. The large 
transition magnetic moment required for a sterile neu- 
trino to decay within a typical detector volume means 
that it can also be an efficient production mechanism. 
This greatly constrains the transition moment with elec- 
tron and especially muon neutrinos. Further study is 
required to investigate whether this production mecha- 
nism could provide exotic but viable explanations of the 
neutrino experiment anomalies. As for the model that 
partially motiviated this study, Ref. [l[ , we find that the 
part of the parameter space with rnjv in the interval be- 
tween 40 and 80 MeV, and the square of mixing angle 
above 10 -3 is in conflict with the RMC constraints. It 
is possible, however, to modify the model and avoid con- 
traints from RMC and magnetic moment induced pro- 
duction by introducing more degrees of freedom beyond 
SM. We find that a radiative decay into another massive 
state, N N'j may provide a remedy. 
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Note added: In a recent comment on our paper [19] . 
Gninenko finds numerically different results from our 
constraints. Our current version includes the effects of 
possible angular anisotropics in the decay of N. We be- 



lieve that the remaining differences between our and Gni- 
nenko's results might be attributed to what appears to 
be only one helicity of N in the final state of OMC in 
Ref. [13, while both should be present. 
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